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ABSTRACT 

Axion-photon mixing has been proposed as an alternative to acceleration as the explanation for 
supcrnovae dimming. We point out that the loss of photons due to this mixing will induce a strong 
asymmetry between the luminosity, cZl(z), and angular-diameter distance, cIa(z), since the latter is 
unaffected by mixing. In a first search for such an asymmetry we introduce a dimensionless mixing 
amplitude A so that A = if no photons are lost and A ~ 1 if axion-photon mixing occurs. The 
best-fit to SNIa and radio galaxy data is A = — 0.31q 4 (95% CL), corresponding to an unphysical, 
negative, mixing length. This same argument limits the attenuation of light from supernovae due to 
dust. We show that future dh and cLa data from SNAP and galaxy surveys such as DEEP2 and KAOS 
will detect or rule out mixing at more than 5a, almost independently of the dark energy dynamics. 
Finally we discuss the constraints from the near maximal polarisation of the gamma-ray burst (GRB) 
GRB021206. Since mixing reduces the polarisation of distant sources, future observations of high 
redshift GRBs may provide orthogonal constraints on axion-photon mixing and related scenarios. 
Subject headings: Cosmology: cosmological parameters 



1. INTRODUCTION 

The reciprocity relation is a wonderfully powerful re- 
sult valid for any metric theory of gravity where pho- 
tons travel on null geodesies, as long as photon number 
is conserved (Etherington 1933; Ellis 1971). It ensures 
that the luminosity distance, di,(z) is exactly the same 
as the angular-diameter distance, cLl(z), up to a factor 
of (1 + z) 2 . In this paper we turn the reciprocity rela- 
tion around and use it to probe alternatives to cosmic 
acceleration. 

The accumulating evidence for recent cosmic accelera- 
tion (Barris et al 2003; Knop et al, 2003) leaves us with 
the familiar coincidence problem - why do we live at such 
a special time? An attractive alternative is that the ac- 
celeration is a mirage and not a real feature of the dy- 
namics of our Universe. 

Although such a non-accelerating cosmology can be 
made reasonably compatible with cosmic microwave 
background (CMB) and large scale structure (LSS) data 
(Blanchard et al. 2003), the dilemma then is to explain, 
without acceleration, the dimming of distant Type la su- 
pernovae (SNIa) and the observed ~ 3cr correlation be- 
tween the CMB and LSS (Boughn & Crittenden 2003; 
Nolta et al. 2003; Scranton et al. 2003; Fosalba et al. 
2003). The latter can perhaps be explained by nega- 
tive spatial curvature, while the dimming of supernovae 
can be explained by axion-photon mixing (Csaki et al 
2002), meaning that the evidence for acceleration is not 
yet overwhelming. 

The basic idea of axion-photon mixing is simple. On 
average, and on large scales relative to the mixing length, 
1/3 of photons in the visible would be lost through con- 
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version to a light axion state, a, in the presence of the 
cosmic magnetic field ~§ . This proceeds through the ax- 
ion interaction term 2^1? ■ ~§. Csaki et al. (2002) argue 
that an axion mass scale of M ~ 4 x 10 n GeU would 
provide a good fit to SNIa luminosity data as a function 
of redshift, (quantified in Erlich & Grojean 2001), with- 
out the need for cosmic acceleration, while still being 
(marginally) consistent with other constraints (see Def- 
fayet et al 2002; Mortsell et al. 2002; Christensson and 
Fairbairn 2003; Mortsell and Goobar, 2003), especially if 
non-flat FLRW models are considered. 

Intriguingly, axion-photon (AP) mixing, with a similar 
mass scale, can explain the existence of super-GZK cos- 
mic rays (Csaki et al. 2003) if the primaries are taken 
to be photons, since they can travel most of the way as 
axions, then oscillate back into photons before reaching 
earth. Axion-photon mixing can therefore provide a si- 
multaneous solution to the super-GZK and coincidence 
problems and is thus worth further detailed study. 

We make four points about the axion-photon mixing 
scenario: 

1. Axion-photon mixing should induce a violation of 
the reciprocity relation and a fundamental dis- 
agreement between the dimensionless coordinate 
distance, y(z), inferred from the luminosity, dziz), 
and angular-diameter distances, dA(z). 

2. The observed dyt(z) should correspond to a de- 
celerating universe if axion-photon mixing is the 
source of supernovae dimming. However, current 
estimates of dA[z), favour an accelerating Universe. 

3. Future data from SNAP and KAOS will allow for 
constraints beyond the ha level. Tests of number 
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counts from the on-going DEEP2 survey will pro- 
vide further tight constraints. 

4. Mixing leads to depletion of the polarisation lev- 
els of extra-galactic sources (Csaki et al. 2002; 
Mortsell and Goobar 2003). The near maximal po- 
larisation seen in the gamma-ray burster (GRB) 
GRB021206 (Coburn and Boggs, 2003) suggests 
that GRB's may provide powerful constraints on 
the mixing scenario when more data is available. 

In this paper we will assume a flat FLRW model consist- 
ing of dust (CIm) and dark energy X (fix = 1 — ^m) 
with a constant equation of state w — px/px- 

2. THE RECIPROCITY RELATION 

One may define several distances in cosmology. The 
luminosity distance, di,(z), estimates distances by com- 
paring the absolute luminosity of an object to its ob- 
served/apparent luminosity. The angular-diameter dis- 
tance, <Ia{z), estimates distances based on how the ap- 
parent linear size of an object changes with redshift. In 
metric theories where photons travel on null geodesies 
and their number is conserved one can show that these 
two distances are fundamentally related by the reci- 
procity relation (Schneider et al, 1992): 



d L (z) = (l + zfd A {z), 



(1) 



When the reciprocity relation holds, the dimensionless 
coordinate distance y(z) can be estimated from either 
<Il(z) or cIa(z) via the relations y(z) = HqcIl/{1 + z) = 
Hq(1a(1 + z), where Hq is the current value of the Hubble 
constant. 

In stark contrast the reciprocity relation is not obeyed 
in the axion-photon mixing scenario, nor indeed in any 
scenario (such as light attenuation due to dust) which 
effectively violates photon number conservation. As a 
result y{z) estimated from dh and dA data should dis- 
agree since d A {z) is unaffected but the luminosity dis- 
tance is modified as d^ — * ^i/y^y^-p where P 7 ^ 7 is 
the probability that a photon will reach earth in a pho- 
ton state and hence be detected. We use the rather good 
approximation 



P -1 + l p -i/i^ 
^7-7-3 + 36 



(2) 



where I is the proper distance of the supernovae. For 
SNIa at cosmological distances the mixing saturates at 
2/3 (Csaki et al. 2002) and hence supernovae should 
appear y/3/2 times further away than they really are, in 
good agreement with the di,{z) predicted by the standard 
best-fit ACDM model with ft A ^ 0.7 and Cl M ~ 0.3. The 
precise value of the decay length /doc depends on the mix- 
ing and the galactic magnetic fields. The preferred value 
of Csaki et al. 2002 is /dec ~ V(2-H"o)- This motivates 
the use of a dimensionless suppression amplitude 



A = (2i7 /dc 



(3) 



which is zero if no mixing occurs and one in the case of 
mixing over cosmological distances. 

3. CONSTRAINTS FROM CURRENT DATA 

If axion-photon mixing is to solve the coincidence prob- 
lem then we should expect the dA(z) data to fit best to 
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Fig. 1. — The residual, Aj/(z), in the dimensionless coordinate 
distance relative to an (Qm,&a) = (1,0) universe for redshift- 
binned SNIa data (di, (2), solid triangles) "corrected" for the loss 
of photons predicted by axion-photon mixing and redshift-binned 
radio galaxy data (cLa(z), circles). The solid curve is Ay(z) for 
a Q.m = 0.22, = 0.78 Universe. If mixing were the origin of 
supernovae dimming we should expect the radio galaxy data to 
coincide with the best-fit y(z) curve of this corrected SNIa data, 
near Ay(z) = 0; however the radio galaxy data lies systematically 
above this curve, favouring an accelerating universe. 



a non-accelerating Universe. This is not the case, how- 
ever. There are (at least) three independent data sets 
for d A {z) that give large best-fit values of consis- 
tent with standard di,{z) best-fits and an accelerating 
Universe. 

Daly & Guerra (2002) and Daly & Djorgovski (2003) 
analysed data for 20 bright FRIIb radio galaxies at red- 
shifts between 0.43 and 1.79 and, assuming a flat uni- 
verse, found fijvf < 0.5 and —2.5 < w < —0.25 at the 
90% confidence level where w is the equation of state 
of the dominant, non-dust component. Another analysis 
(Jackson, 2003; Jackson and Dodgson, 1997) of ultra- 
compact radio sources (Gurvitis 1994; Lima and Alcaniz 
2002) at z > 0.5 found that the best-fit flat ACDM model 
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Searches for co-moving standard rulers via peaks in 
the two-point correlation function of quasars have also 
been undertaken. Mamon and Roukema (2002), using 
a subset of the 2QZ 2df quasar survey, estimated Qa — 
0.65 ± 0.35. Assuming a flat universe they constrain the 
equation of state of the non-dust matter to w < —0.35 
at 2cr. This approach has been extended recently by 
Outram et al (2003) using the full 2QZ survey, allowing 
for even stronger results. Assuming a flat FLRW model, 
they find Q\ — 0.71*®' ^ and exclude an Q\ = Universe 
at over 95% confidence. 

In summary, all current estimates of g?a(-z) favour an 
accelerating Universe, and since they arc unaffected by 
axion-photon mixing, disfavour it as the explanation for 
the majority of SN-Ia dimming. This point is made 
visually in figure Q which shows the binned valus of 
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Fig. 2. — The 2d likelihood plot for the equation of state w 
and Qm for the combined SNIa & RG data sets, with 1 and 2<r 
contours shown. A flat FRW universe was assumed. The combined 
data clearly favour an accelerating low-density universe. 



the dimensionless coordinate distance, y{z), of the data 
sets used in our analysis. The SN-Ia data has been cor- 
rected for mixing, and indeed clusters around the non- 
accelerating EdS model (17m, ^a) = (1,0). We used the 
combined data sets of Tonry et al (2003) and Barris et 
al (2003) to which we added the new supernovae from 
Knop et al (2003). However, the radio galaxy data clearly 
prefers an accelerating model. This data is a combina- 
tion of Daly and Djorgovski (2003), Jackson (2003) and 
Gurvits (1994). A more detailed description of the data 
sets can be found in Bassett and Kunz (2003). 

In figure J2J) we show the 2d likelihood plot for Q,m and 
w for the combined data sets assuming a flat Universe. 
The overall best-fit with axion-photon mixing is VLm = 
0.24, w = —1.1 and A = —0.3. The one-dimensional 
(marginalised) 95% confidence limits on the parameters 
are 0.15 < ft M < 0.33, -1.6 < w < -0.6 and -0.7 < 
A < 0.3. The values preferred by axion-photon mixing 
(w = — 1/3 and A ~ 1) are ruled out at well over 3er. 

Negative values of A correspond to the appearance of 
photons, instead of absorption. This is not impossible 
in the axion-photon scenario, e.g. if SN-Ia produce large 
numbers of axions which then become photons on the 
way to earth. Still, it may be argued that this region 
of the parameter space is unphysical and should be ex- 
cluded. In this case the 95% upper limit (one sided) on A 
is 0.6, and A = 1 lies at 3a. This is no longer sufficient to 
reliably rule out axion-photon mixing. But the limits on 
the equation of state remain rather strong, with a 95% 
confidence interval of —1.1 < w < —0.5 and w = —1/3 
remains ruled out at over 3a. This renders the scenario 
quite unattractive, as the mixing does not alleviate the 
coincidence problem asssocitated with cosmic accelera- 
tion. 

4. CONSTRAINTS FROM FUTURE DATA 

Future estimates of (Il(z) from the SNAP satellite 1 
and (1a{z) from number counts from the DEEP2 survey 

1 http://snap.lbl.gov 
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Fig. 3. — A — w likelihood plot for current SNIa and RG 
data which peaks at (A,ui) = (—0.3,-1.1), favouring no mixing 
(A = 0) and acceleration. Axion-photon mixing corresponds to 
A fa 1. The likelihood contours for future data, based on a fiducial 
non-accelerating model with mixing: (A, w) = (1,-1/3), are also 
shown, at the top right. 




Fig. 4. — The Id likelihood plots for A and w for current data 
(solid line) and hypothetical future data (dashed curve). A = 
corresponds to no photon loss, A = 1 corresponds to axion-photon 
mixing. Current data gives A = —0.3 ± 0.25 (lcr). The fiducial 
model for the future data is flat with (A. w) = (1, —1/3) and Hm = 
0.3. 



and from baryon oscillations from the KAOS survey 2 
will allow estimates of Vl{z) and ua(z) at the level of a 
few percent (Aldering et al, 2002; Linder 2003; Blake & 
Glazcbrook 2003; Seo & Eisenstein 2003). To investigate 
the power of future experiments we used the errors given 
in the dark energy science case of the KAOS purple book 
for both the KAOS and SNAP experiments. The central 
values of the data are chosen to match a model with 
mixing with an underlying flat FLRW cosmology with 
= 0.3 and w = —1/3 which corresponds to the best- 
fit of Csaki et al. (2002). 

Assuming the auxiliary cosmic parameters (e.g. fljf) 
are well-known from other methods by then, we halved 
the current best estimates and assumed Qjvf = 0.3 ±0.02 
as a prior. Although this is not required (and indeed our 
analysis with current data does not make this assump- 
tion), it helps to reduce the error on w significantly, and 
it certainly is sensible 3 . With these assumptions we con- 

2 http://www.noao.edu/kaos/ 

3 When comparing our errors on the equation of state with those 
of the KAOS purple book, one should note that they fit simultanc- 




4 



elude that we will be able to detect or rule out the mixing 
scenario at over 5a after marginalising over w. The esti- 
mated error bar on A is less than 0.1 (and is degraded to 
about 0.13 if no constraints on f_M are added). The 2d 
likelihood in the A — w plane is shown in the upper right 
hand corner of Figure J3J • 

Although we have assumed a constant w here, the 
beauty of having both cLl and (La information is that it al- 
lows us to separate the issue of mixing from the dynamics 
of the dark energy. At each redshift, axion- photon mix- 
ing should lead to fundamentally inconsistent values for 
y(z) derived from <_£ and d,A respectively. 

An estimate of cosmic parameters unbiased by mixing 
will be available from the Sloan Survey (Matsubara & 
Szalay 2003) while a further test of mixing is provided 
by number counts versus redshift, dN/dz, which depends 
on dA(z). Since the volume of space as a function of 
redshift is very sensitive to A number counts is a good 
test of acceleration. Generally, the number of objects in 
the range of affine parameter values [y,y + Ay], is (e.g. 
Ellis 1971, Ribeiro k, Stoeger 2003) 



dN 



l(1 + z) n(y) dy dfl 



(4) 



where n(y) is the number density of objects and dfl is 
the differential solid angle at the observer. Axion-photon 
mixing alters galaxy number counts by reducing the ap- 
parent luminosity of objects at high redshift, at least in 
the visible range. Since high redshift objects appear dim- 
mer, the selection function ip is altered and faint galaxies 
will be lost. Therefore, there should be a deficit of ob- 
jects relative to the case of no axion-photon mixing. 

If we therefore compare a standard ACDM model 
against a non-accelerating model with axion-photon mix- 
ing the difference in number counts at z > 1 is significant. 
One may consider variants of this basic idea such as the 
dV/dzdQ. test which, applied to the DEEP2 galaxy sur- 
vey of ~ 50/c galaxies with redshifts 0.7 < z < 1.4, should 
allow an estimate of w today (unbiased by axion-photon 
mixing) to ~ 10% (Newman & Davis, 2000). 

The mixing mechanism may be constrained in yet an- 
other manner however. Observations of the polarisa- 
tion of light from gamma-ray burst (GRB) GRB021206 
(Coburn and Boggs, 2003) have found linear polarisa- 
tion levels of II = 0.80 ± 0.2, centered very near the 
maximum allowed by Compton scattering which strongly 

ously several parameters more, and that our fiducial model has a 
higher value of w = —1/3, which helps to reduce the errors further. 



supports synchrotron radiation as the source of at least 
some GRB's. If GRB021206 is at a redshift z > 0.1 and 
Compton scattering is the source of the linear polarisa- 
tion, then the near maximal value of II observed on earth 
leaves little room for depletion due to mixing. However, 
as pointed out in Csaki et al. (2002), mixing is intrinsi- 
cally inhomogeneous. It is possible to have certain lines 
of sight that experience essentially no mixing at all, de- 
pending on the magnetic field traversed. Hence, unless 
there is a high-z SNIa in the same narrow field of view 
as the GRB, a single event alone cannot rule out the 
mixing scenario. Further, the linear polarisation of the 
GRB may not be due to Compton scattering (Lazzati 
et al, 2003), in which case there might still be room for 
axion-photon mixing. 

5. CONCLUSIONS 

The dimming of distant supernovae (SNIa) remains the 
most direct evidence for cosmic acceleration. Neverthe- 
less alternative explanations exist, such as axion-photon 
mixing in which roughly one third of all photons from 
distant SNIa are lost into axion states. We have pointed 
out that such mixing will not affect the angular-diameter 
distance <-a(z) and hence will cause a fundamental asym- 
metry between measurements of the luminosity distance, 
di,(z), and d^(z) that can be searched for. 

In a first search for such asymmetry we have under- 
taken a joint analysis of high-redshift SNIa (c?l(z)) and 
radio galaxy data (<Ia{z)). The results do not favour the 
loss of photons and hence disfavour mixing. Future data 
will improve the limits, and be able to test very generally 
for an asymmetry between <1a(z) and dt{z)- Number 
counts versus redshift are a promising test while esti- 
mates of c?l(z) from SNAP and cU(z) from a large 2nd 
generation galaxy survey such as KAOS will allow axion- 
photon mixing to be detected or ruled out at more than 
5<7, almost independently of the dynamics of the dark 
energy, showing the power in constraining non-standard 
physics implicit in combining and g?a(z) data. 
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